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BOPACE THEORETICAL MANUAL - ADDENDUM 


2.6 IMPROVED ALGORITHM TOR INELASTIC CALCULATIONS 


Suwnary of Basic Concept s - The iterative residual -force procedure Is 
often employed with an incremental solution for inelastic (plasticity 
and creep) problems, in order to avoid accumulated error. Each iteration 
in the residual -force procedure involves the following two stages. 

1) Equilibrium and Compatibility: Given the current residuals (unbalanced 

forces or stresses), the equilibrium and compatibility equations 

are applied in order to predict an improved configuration (of 
displacements and strains). 

2) eparation of Strains: Given the current strains, some algorithm 

based on the inelastic material theory is applied in order to 
separate the strains into their elastic, plastic and creep portions, 
and thus provide fie resulting stresses. 


) 
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When this procedure has converged to the correct result, the following 
conditions will be met. 

1) Forces In equilibrium 

2) Displacements compatible 

3) Plastic strain increment satisfies normality rule 

4) Size of yield surface consistent with deformation history ^ 

5) Translation of yield surface consistent with deformation history 

The overall BOPACE solution technique based on the residual -force 
procedure is summarized in Section 4. The purpose of the present section 
is to discuss the details of a new algorithm which has been developed 
and Incorporated into BOPACE, for improving the convergence and accuracy 
of the inelastic stress-strain calculations. This algorithm defines the 
implementation of stage 2 (separation of sv.rains) in the residual-force 
iterative procedure. 

B ackgroun d - The theory already presented in Sections 2.1 through 2.5 
may be employed for both stages of the iterative procedure, and in fact 
equations of the type 2.5-5 were used for all stress-strain calculations 
in the initial version of BOPACE. Convergence difficulties resulted 
from the use of this approach in stage 2, however, when the incremental 
inelastic strains were large relative to the cumulative elastic strains. 
These difficult : es were substantially eliminated by properly controlling 
the direction defined for the incremental inelastic strains. (The 



reason for the difficulties and the method of control were presented in 
Reference A-l). A.iother quite different approach is based on a "strain- 
space" concept, and was presented by Barsoum in Reference A-2 wi th the 
claim of a significant impro ement in efficiency. Ti.it approach has 
therefore also been developed and evaluated for BOPACE usage. Because 
the method as presented in Reference A-2 assumes kinematic hardening 
only, it was extended to incorporate the combined isotropic and kinematic 
h 'dt.-'ing provided by BOPACE. In addition, some further techniques for 
accelerating convergence were identified and incorporated into the 
method. The resulting modified version of BOPACE has shown promising 
results. For example, a small test problem involving 15 load Increments 
and solved using BOPACE solution option 5, required 25" fewer iterations 
and 50% fewer stiffness matrix updates than when run on the previous 
version. Although farther comparisons are needed to prove the effective- 
ness of the new algorithm, it is expected to supersede the previous 
BOPACE capability. 

Basic Definitions a nd Compar ison o f Algorithm s - The new inelastic 
algorithm involves calculations in the "deviatoric strain space", rather 
fhan the more conventional "deviatoric stress space" used 1 r. previous 
BOi ACE programs. For the sake of clarity, the previously used stress- 
space algorithm wi.l again be summarized here, and the elastic-plastic 
quantities used in the new strain-space algorithm will be defined and 
compared with previous quantities. 




As described in Section 2.3, the definition of a plasticity theory 
requires assumptions for three basic constituents: a yield surface, a 

flow ru’e, and a hardening assumption. BOPACE development is based on 
the Mises yield surface, and this surface is represented by a hypercircle 
in 9-dimensional deviatoric stress space, as shown in Figure 2.6-la. 

The surface is defined by thj equation 
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w^re s is the deviatoric stress, s = s - a is the relative deviatoric 
stress and defines the isotropic hardening, i is the surface translation 
and defines the kinematic hardening, while s° is a reference value of s 
and must be known as a function of plastic deformation (e.g. frou a 
uniaxial test). Point A in Figure 2.6 la is the origin of the ceviatoric 
stress space, point B is the current center of the yield surface, and 
point C represents the current state of deviatoric stress. A stress 
point on the surface corresponds to a plastic ;t. te. According to the 
Prandtl -Reuss flow rule, the direction of tht emental plastic strain, 
Ae P , is normal to the yield surface at the current deviatoric stress 
ite, s. A solia circle ( • ) in Figure 2.6-1 denotes a point which 
remains fixed throughout the increment, while an open circle ( 0 ) denotes 
a point which moves during the increment. In order to achieve greater 
accuracy and allow larger load increments, BOPACE evaluates moving points 
such as B and C at the midpoint of the plastic increment. Additional 
details of the BOPACE stress-space algorithm are discussed in Section 23. 
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For the new strain-space algorithm, the three basic constituents of the 
plasticity theory remain unchanged, and direct use is made of the stress- 
space theory and nomenclature. However, we now work with a yield surface 
and associated quantities in strain-space. Thus we compute the deviatoric 
elastic strain, e P , in terms of the deviatoric stress, s, by 

e* - s t /G (2.6-?) 

f 

where G = E/(l+v) is a tensorial shear modulus. Similarly we define a 
"strain center", n, in terms of the stress center, ., by 

6 i -<yG (2.6-3) 

Then the relative deviatoric strain, e. Is defined by 

e 1 B e* - ■ (s i - a < ) /G - Sj/G (2.6-4) 

The geometrical interpretation of the new algorithm Involving these 
o^antities is provided by a sketch in 9-dimensional deviatoric strain- 
space, shown in Figure 2.6-lb. There point 0 is the origin, defining 
/ *> initial undeformed (zero strain) state. Subsequent deformation is 
caused by a series of load increments, resulting in elastic and plastic 
strains. A superscript 0 is used to denote the value of a quantity at 
the beginning of the load increment. Thus, point A defines the cumulative 
plastic strain, e^'\ which exists at the beginning of the current increment. 
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(Because of the plastic incompressibility assumption, the plastic 
strains themselves are deviatoric strains). All other pcints in Figure 

2.6- lb refer to locations at some time during the current increment. In 
particular, we will be mainly concerned with the location of these 
points at a defined reference time. This reference time may be taken at 
the end of the increment, following the approach of Barsoum [14], or 
greater accuracy may be obtained at the expense of some additional 
variable storage by taking the reference time at the midpoint of the 
plastic Increment, as Is done in the new BOPACE algorithm. Point 0 
defines the total cumulative deviatoric strain, e, during the increment. 
The circle is associated with the Mises yield surface, but is a hyper- 
circle in the deviatoric strain space. A strain point withir the 
surface corresponds to an elastic state, while a strain point outside 
the surface corresponds to a plastic state. The size of this circle is 
defined by its radius e^e^ 3 s^/G), whereas the Mises stress-space 
surface has radius s. The center of the circle is at point B (B. 3 

e P ^ + + a j/G)« whereas the center of the Mises stress-space 

surface has components a • . During plastic deformation, the strain-space 
surface may undergo both expansion (due to isotropic hardening), and 
translation (due to kinematic hardening). The cumula;ive deviatoric 
plastic strain, e e , is defined by the vector AC (e*: = s^/G). From these 
comparisons it should be apparent that the basic quantities in Figures 

2.6- la and b, respectively, can be made to coincide, if points A 
are superimposed and all dimensions in 2.6-lb are divided by the 
factor G. The incremental plastic strain, At P , is defined by the 
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vector CD. It is normal to the circle because of the Prandtl -Reuss 
flow rule, and is therefore col inear with the radius e to point C. 

The vector BD * e ♦ ^ is denoted by E‘. The symbols e', e^, Ac* 5 

and E' are consistent with their usane in Reference A-2. 

Computation Procedure - We now define the new strain-space algorithm 
for implementing stage 2 r f the res idua 1 -force iterative procedure. 

The problem which must be solved can be stated in terms of the various 
strain vectors. At the beginning of the increment, we have known 
values for (which remains constant during the increment), and for 
B, e e , and e. These have been determined such that they are all 
consistent, i.e., such that the appropriate veci -s meet at single 
points A, B and C. The current estimate for the value of e' at thp 
reference time is also known from stage 1 of the iterative procedure. 

We must then determine values for tf, e , e and Ac k at che reference 
time, consistent with the convergence requirements. Stated somewhat 
differently, we are given the locations of points A and D at the reference 
time, and the locations of points B and C at the beginnlna of the increment. 
We must then compute the locations of B and C at the reference time, 
consistent with the convergence requirements. 

The basic steps of the stage 2 algorithm are summarized by the following. 

1) Given values at beginning of increment for: 

0 

a » 


stress center 



0 

s * relative deviatoric stress 

of) 

e • elastic strains 

Given Ac ■ total (elastic ♦ plastic) strain Increment from staqe 1. 
Compute values a f eference time, based on estimated Incremental 
deformation, for: 

Aa « kinematic hardeninq Increment 

AS ■ isotropic hardeninq Increment 

eO * 

Compute: ♦ Ac^ ■ Initial elastic strain ♦ total strain 

increment 

I 

e. « correspondinq deviatoric value 
Compute: 8 j * (<«^ + Aa^)/G 

e. * ( ► as^/G 

Compute E! 3 ej - 6^ 

Compute \ 8 ( | E ' | - |e|)/|E'| * plastic proportionality constant 

Compute AeP = aE! = Incremental plastic strain at reference time. 

Adjust AcP «- At.P times (ratio of total to reference time incre- 
ment), to obtain total plastic strain increment. Compute 
Ac® * Acj - AeP * incremental elastic strain. 

Compute end of increment values for: 

c® * r.®^ + At® •* cumulative elastic strain 
c. * e® * cumulative stress 

Use o to compute residual forces and error norm, and return to 
staqe 1 if converqence has not been achieved. 
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The strain-space algorithm presented above corresponds to that given by 
Barsour [A-2J except that here a combined isotropic and kinematic hardening 
is provided and a reference (midpoint) time calculation of the incremental 
variables is used to improve accuracy. As noted by Barsoum, greater 
consistency and better convergence are obtained by utilizing an algorithm 
In strain space rather than in stress space. This Is because the stress- 
space calculation fixes the Ai P vector along the direction of the current 
s vector, rather than simultaneously fixing the directions of s and 
Ac P consistent with the given total strain Increment Ac. The stress- 
space Iteration can cause large oscillations In the location of point C, 
resulting in divergence if At P is large relative to the cumulative 
elastic strain. 

Although a strain- r pace algorithm eliminates most of the inconsistencies 
and tendencies toward divergence, it should be noted that an incon- 
sistency still exists in the plastic hardening quantities. This Is 

A 

because Aa and As are based on the estimated incremen . of plastic defor- 
mation, which will not in general be consistent with the actual deforma- 
tion. Thus if another iteration were performed using the same value for 
the total strain increment Ac, different results would bo obtained due 
to change in 6 and e. 
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This difficulty Is eliminated in the present approach by properly 
modifying the calculation of a in step 7. For this calculation, we 
use the parameters c and r associated with kinematic and isotropic 
hardening, respectively, >•’. the expressions 


**1 " ! c At ? 

ASi - | r Aef 
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We then see that 

E i " e i ‘ 8 i " e i - (0° ♦ 

• e i - 0? “ Acyc - ej - 00 _ 2 c Atf p /G ( 2 . 6 - 6 ) 


Replacing Ae? in this equation by AEj , we may solve for E!: 


E 1 " te i ’ B?)/(l + f A c /G) (2.6-7a) 
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In a similar manner we may obtain 


e 1 ■ e° ♦ v A r tJ/G 

The plastic proportionality constant, as already defined, is 

A - (|e'| - Je | )/ |C ' 1 (2.6f7c) 

It is apparent from Equations 2.6-7 that the expression for a Is non- 
linearly dependent upon a Itself, and this is the reason why a con- 
sistent A is not solved for directly. An accurate value for x, however, 
can easily be obtained by a 'linear Intersection method." In this 
method we take the approximate value of x from step 7, and substitute 
into the Equations 2.6-7 to obtain a new computed value x^. We then 
assume a value of A + ax, where ax is a small change (perhaps .01 x), and 
again substitute into Equations 2.6-7 to compute another value x^. 

The two pairs of assumed and computed x values are plotted in Figue 
2.6-2. The correct value for x lies on the 45-degree line (since there 
the assumed and computed values would be equal), at the intersection of 
♦his line with the line connecting the two plotted points. This cor- 
rected value of X is obtained by the following adjustment of x from step 
7. 
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\ ♦ X ♦ AAU - * c0 )/( x c i " A cO _AA ) 


The Incorporation of this adjustment into the strain-space algorithm 
provides consistent values for all quantities in stage 2 of the iterativ 
process, ano results in improved convergence. 



ASSUMED VALUE (A) 


Figure 2.6-2. Linear- Intei section Calculation fc- \ 
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Extensions and Refinement', to the basic Algor ith m - The strain-space 
alqorithm as presented here Is employed in BOPACE for plastic analysis. 

In addition, the ROPACE alqorithm treats creep strains In a manner 
similar to that *or the plastic strains, for cases where the material 
Is elastic at the beqinninq of an increment and then reaches the plastic 
yield point at some intermediate f ie durinq the increment, qreater 
accuracy is obtained by dividing the calculations into two parts. In 
such cases the initial creep is taken in the direction of the initial 
deviatoric stress, and creep which occurs after the yield point is taken 
in the same direction as the plastic strain increment. Other refinements, 
such as temperature dependent elastic-plastic-creep and qeneralized load 
reversal, are treated as discussed in Section 2. 
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